Quasinormal modes and second order thermodynamic phase 
transition for Reissner-Nordstrom black hole 



Jiliang Jingj and Qiyuan Pan 

Institute of Physics and Department of Physics, 
Hunan Normal University, Changsha, 
Hunan 410081, P. R. China 
and 

Key Laboratory of Low Dimensional Quantum Structures 
and Quantum Control (Hunan Normal University), 
Ministry of Education, P. R. China 

Abstract 

The relation between the quasinormal modes (QNMs) and the second order thermodynamic phase 
transition (SOTPT) for the Reissner-Nordstrom (RN) black hole is studied. It is shown that the 
quasinormal frequencies of the RN black hole start to get a spiral-like shape in the complex u plane 
and both the real and imaginary parts become the oscillatory functions of the charge if the real part 
of the quasinormal frequencies arrives at its maximum at the second order phase transition point of 
Davies for given overtone number and angular quantum number. That is to say, we can find out the 
SOTPT point from the QNMs of the RN black hole. The fact shows that the quasinormal frequencies 
carry the thermodynamical information of the RN black hole. 
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The QNMs of a black hole are defined as proper solutions of the perturbation equations 
belonging to certain complex characteristic frequencies which satisfy the boundary conditions 
appropriate for purely ingoing waves at the event horizon and purely outgoing waves at infinity 
1|. They are entirely fixed by the structure of the background spacetime and irrelevant of 
the initial perturbations . Thus, it is generally believed that QNMs carry the footprint 
to directly identify the existence of a black hole. Meanwhile, the study of QNMs may lead 
to a deeper understanding of the thermodynamical properties of black holes in loop quantum 



gravity 



3| 4j, as well as the QNMs of anti-de Sitter black holes have a direct interpretation in 



terms of the dual conformal field theory {5], f], H] . 

On the other hand, one of important characteristics of a black hole is its thermodynamical 
properties. It is well known that the heat capacity of the Schwarzschild black hole is always 
negative and so the black hole is thermodynamical unstable. But for the RN black hole, 
the heat capacity is negative in some parameter region and positive in other region. Davies 
pointed out that the phase transition appears in black ho 
takes place at the point where the heat capacity diverges 

Because the QNMs of a black hole are entirely fixed by the structure of the background 
spacetime and the SOTPT is only related to the parameters of the black hole, an interesting 
question is whether there are some relations between them. The aim of this paper is to study 
this question for the RN black hole, and we find that the quasinormal frequencies do carry the 
thermodynamical information of the black hole. 

Assuming that the azimuthal and time de pen dence of the fields will be the form e - l M- m </?) 



e thermodynamics and the SOTPT 

3, H, H ■ 



and using the Newman- Penrose formalism 



we can obtain the separated equations for 
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massless scalar, Dirac, and Rarita-Schwinger (RS) perturbations around a RN black hole 

[A£>i_ s £>J + 2(2s - l)iur - (A + 2s)}A s R s = 0, 

\C\_ a C 8 + (A + 2s)}S s = 0, (s=0, 1/2, 3/2) (1) 



[AV\ +S V + 2(2s + l)iur - X)R S = 0, 

(C 1+s Cl s + X)S s = 0, (s=0, -1/2, -3/2) (2) 
where A = r 2 — 2Mr + Q 2 , M and Q represent the mass and charge of the black hole, and A 
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is the angular separation constant [la. lid. Il7| 

(l-s){l + s + l), l = \s\,\s\ + 1 r~ , 

\ & ) 

(j - + s + l), j = \s\, \s\ + 1, • - • , 
where I and j are the quantum number characterizing the angular distribution for the boson and 
fermion perturbations respectively. Introducing an usual tortoise coordinate dr* = r 2 /Adr and 
resolving the equation in the form R s = A~ s / 2 \I/ s /r, we can rewrite the radial wave equations 
in Eqs. © and © as 

d 2 ^ t 
dr\ 

with 



+ [lu 2 -V]V s = 0, (4) 



V = lS .r 2 ^ + (^ + A ) A +(ff) 2 + 

dr r 4 r 4 r 3 dr r 2 

The boundary conditions on wave function at the horizon and infinity can be expressed 

as 

(r - r+) 2 2K + r -> r+, 

^ . (6) 

where k± = (r+ — r_)/(2r 2 t ) is the surface gravity on the horizons r±. A solution to Eq. (|4]) 
that has the desired behavior at the boundary can be written as 

s iu> 1 S i o„' I 

t / \ ~~ o ~~ o / , 1 ^ -\-Z LUJ-\--x 

oo ✓ \ m 

xe-<"->E« m (^) • (7) 

m=0 ^ ' 

If we take r + + r_ = 1, the sequence of the expansion coefficients {a m : m = 1,2,....} is 
determined by a three-term recurrence relation staring with do = 1: 

a ai + p a = 0, 

a m a m+1 + (3 m a m + 7 m a m _i = 0, m=l,2,.... (8) 

The recurrence coefficient a m , j3 m and 7 m are given in terms of m and the black hole parameters 
by 

a m = m 2 + (C + l)m + C , 
/3 m = -2m 2 + (d + 2)m + C 3 , 

7m = m 2 + (C 2 - 3)m + C 4 - C 2 + 2, (9) 



and the intermediate constants C m are defined by 



Cq = 1 — s — ioo — iuB, 
d = -4 + 2iu(2 + b) + 2iujB, 
C2 = s + 3 — 3iuj — iuB, 
C 3 = uj 2 (A + 2b - 4r+r_) - s - 1 + (2 + 
-A 2 + (2w + 

C 4 = s + 1 + 2 ( iuo - s - - J za; - (2cj + i)wB, (10) 
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where -B = (r 2 ^ + r 2 _)/{r + — r_). The series in ((7]) converges and the r = +00 boundary 
condition (jSJ) is satisfied if, for a given s and A, the frequency u is a root of the continued 
fraction equation 

Otm-llm <y m -2lm-l ^071 



(m = l,2..X (11) 



This leads to a simple method to find quasinormal frequencies of the RN black hole — defining 
a function which returns the value of the continued fraction for an initial guess at the frequency, 
and then use a root finding routine to find the zeros of this function in the complex u plane. 
The frequency for which happens is a quasinormal frequency [3, \\\ . 

The Figs. EH describe the QNMs for the scalar (s=0), Dirac (s=-l/2), and RS fields (s=-3/2) 
obtained by the continued fraction method respectively. In Fig. [T] (or 121125 f° ur P ane ls show 
trajectories in the complex uj plane of the scalar (or Dirac, RS) quasinormal frequencies of the 
RN black hole with different angular quantum number and overtone number. The other panels 
draw the imaginary part Im{uo) and real part Re(uj) of the quasinormal frequencies versus the 
charge Q. These figures tell us that in the complex uj plane the quasinormal frequencies will get 
a spiral-like shape as the charge Q increases to its extremal value when the overtone number 
equals to or exceeds a critical value n c for a fixed angular quantum number (say, n c = 1 with 
I = and n c = 4 with I = 1 for the scalar field, n c = 2 with j = 1/2 and n c = 5 with j = 3/2 
for the Dirac field, and n c = 5 with j = 3/2 and n c = 9 with j = 5/2 for the RS field.), and at 
same time both the real and imaginary parts become the oscillatory functions of the charge. 
The critical value of the overtone number n c increases as the angular quantum number I (or 
j) increases for a given perturbation, and it also increases as the absolute value of spins \s\ 
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FIG. 1: Left four panels show trajectories in the complex uj plane of the scalar quasinormal frequencies 
of the RN black hole for 1 = 0, n = 1, 2 and 1 = 1, n = 3, 4. The other panels draw the imaginary 
part Im(uj) and real part Re(uj) of the quasinormal frequencies versus the charge Q. 

increases for the same level of the angular quantum number (here the least level is defined as 
I = 0, j = 1/2, and j = 3/2 for the scalar, Dirac, and RS fields respectively). 

On the other hand, we know that the heat capacity Cn of the RN black hole is given by 

m 

TS 3 M 



a 



Q 



(12) 



\dTj Q nQ 4 /4-S 3 T 2 ' 

where T is the temperature and S is the entropy of the black hole. It is obvious that the singular 
point of the heat capacity (SPHC) occurs at ttQ 4 /4-5 3 T 2 = 0, i.e., Q sp = V3M/2 w 0.4330127. 
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FIG. 2: Left four panels show trajectories in the complex uj plane of the Dirac quasinormal frequencies 
for j = 1/2, n = 1, 2 and j = 3/2, n = 4, 5. The other panels draw the imaginary part Im(u) and 
real part Re(uj) of the quasinormal frequencies versus the charge Q. 

Fig. H] shows that the heat capacity is negative in the region Q < Q sp and positive in the region 
Q > Qsp- Davies , Q showed that the SPHC is just the SOTPT point of the black hole 
thermodynamics. 

To find the relation between the QNMs and the SOTPT point, we present critical point 
Q cp , which is the value of the charge at which the real part of the first oscillatory quasinormal 
frequencies arrives at its maximum, with the different critical overtone number n c for the scalar, 
Dirac and RS fields in Table [H From the table we find that the SPHC Q sp is in good agreement 
with the critical point Q cp because the difference between Q sp and Q cp is less than 2.55%. 
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FIG. 3: Left four panels show trajectories in the complex u plane of the RS quasinormal frequencies 
for j = 3/2, n = 4, 5 and j = 5/2, n = 8, 9. The other panels draw the imaginary part Im{uS) and 
real part Re(uj) of the quasinormal frequencies versus the charge Q. 



Besides, for the Schwarzschild black hole there is no SPHC because its heat capacity is always 
negative and there is no critical point because its QNMs never get a spiral-like. 

To compare them farther more, we take K = as the slope of the QNMs. For Q < Q cp , 
Fig. [5] shows that the QNMs have a positive slope and hence large values of ujr correspond to 
large values of uj. However, for Q > Q cp , the slope is negative and hence large values of ojr 
correspond to small values of ujj. At Q = Q cp , the slope becomes infinity. The transition from 
negative to positive value of K occurs via a infinite discontinuity, characteristic of a second 
order phase transition. The trajectories of the slope of the QNMs take the same form as that 
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FIG. 4: The panel shows trajectory of — Cq versus Q near the singular point of the heat capacity 
Q sp = 0.433. 



TABLE I: The values of the charge for Q 



cp ■ 





Scalar field (s= 0) 


(I, n c ) 


(0, 1) 


(1, 4) 


Qcp 


0.438 


0.432 


|AQ| 
Qcp 


1.15% 


0.23% 




Dirac field (s= -1/2) 


(j, n c ) 


(1/2, 2) 


(3/2, 5) 


Qcp 


0.442 


0.436 


|AQ| 
Qcp 


2.08% 


0.69% 




Rarita-Schwinger field (s= -3/2) 


(j, n c ) 


(3/2, 5) 


(5/2, 9) 


Qcp 


0.444 


0.430 


|AQ| 
Qcp 


2.55% 


0.69% 



of — Cq versus Q. 

We all know that there are two characteristic parameters for any perturbation of a black 
hole background: the oscillation time scale r R = 1/ur and the damping time scale 77 = 1/|^/|- 
Near the critical point, although the increment of 77 increases monotonously as Q increases, the 
increment of tr decreases as Q increases for Q < Q cp and it dose not change at all at Q = Q cp , 
but it increases as Q increases for Q > Q cp . It is curious that the change of the increment of 
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FIG. 5: The panels show trajectories of the slope K = dui/diVR versus Q near the critical point Q cp 
for the first oscillatory QNMs of the scalar, Dirac, and RS fields. 

tr presented here is similar to that of the temperature of the black hole when it takes the same 
quality of heat. If we take a slope of the time scale as K T = drj/drn, the trajectory of the slope 
takes the same form as that of Cq versus Q because K T = —K^ujr/uj^ 2 . Above discussions 
show us that the critical point of the QNMs may be associated with the SOTPT point. 

In summary, we study the relation between the QNMs and the SOTPT for the RN black 
hole and find the following results: if the real part of the quasinormal frequencies arrives at 
its maximum at the SOTPT point for given overtone number and angular quantum number, 
the QNMs will start to get a spiral-like shape in the complex uj plane, and both the real and 
imaginary parts will become the oscillatory functions of the charge. The QNMs will (not) take 
a spiral-like shape if the angular quantum number or the overtone number larger (less) than 
this given value. If a black hole does not possess the SOTPT point, its QNMs never take a 
spiral-like in the complex uj plane. Besides the fact that the critical point Q cp of the QNMs 
is in good agreement with the SOTPT point Q sp , the transition from negative to positive 
value of the slope K (K T ) of the QNMs (time scale) occurs via a infinite discontinuity, and 
the trajectories of K and K T take the same form as that of the heat capacity. These facts 



9 



show that the critical point of the QNMs may be associated with the SOTPT point and the 
quasinormal frequencies carry the thermodynamical information of the RN black hole. 
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